clc;
clear;

syms x y z t
u1 = -2*(1-cos(2*pi*x))*sin(2*pi*y)*sin(2*pi*z)*cos(t);
u2 = sin(2*pi*x)*(1-cos(2*pi*y))*sin(2*pi*z)*cos(t);
u3 = sin(2*pi*x)*sin(2*pi*y)*(1-cos(2*pi*z))*cos(t);

B1 = 2*sin(pi*x)*cos(pi*y)*cos(pi*z)*cos(t);
B2 = -cos(pi*x)*sin(pi*y)*cos(pi*z)*cos(t);
B3 = -cos(pi*x)*cos(pi*y)*sin(pi*z)*cos(t);

w1 = (1-cos(2*pi*x))*sin(2*pi*y)*sin(2*pi*z)*cos(t);
w2 = sin(2*pi*x)*(1-cos(2*pi*y))*sin(2*pi*z)*cos(t);
w3 = sin(2*pi*x)*sin(2*pi*y)*(1-cos(2*pi*z))*cos(t);

p = (sin(4*pi*x)+sin(4*pi*y)+sin(4*pi*z))*cos(t);
